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ABSTRACT 

We report new limits on the Gunn-Peterson effect at a redshift near 5 us- 
ing spectroscopic observations of the z — 5 Sloan Digital Sky Survey quasar, 
J033829.31+002156.3, made with the LRIS and HIRES spectrographs on the 
Keck telescopes. Lower resolution spectrophotometrically cahbrated observations 
made with LRIS over the wavelength region AA4500 — 9600 A were used to obtain 
a continuum shape and to flux calibrate much higher resolution {R — 36, 000) 
observations made with HIRES. The LRIS data show an Oke Da index of 0.75. 
Portions of the HIRES spectrum return to near the extrapolated continuum level. 
Including both statistical and systematic errors we place an upper limit of r = 0.1 
on the regions of minimum opacity. We argue that, even if this opacity arises in 
underdense regions of the intergalactic gas, we require a high value of the meta- 
galactic ionizing flux at these redshifts {J^, ^ 4 x 10"^"^ ergs cm~^ s~^ Hz~^ sr~^ 
at 2; ~ 4.72) to produce a solution which is consistent with even minimum nucle- 
osynthesis estimates of the baryon density. We also report the presence of a Mg II 
absorption system of extremely high equivalent width (WA,rest(2796) = L73 A) 
at ^ = 2.304. 

Subject headings: cosmology: observations — early universe — intergalactic 
medium — quasars: absorption lines 
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1. Introduction 

Soon after the discovery of cosmologically distant objects it was realised (Field 1959) 
that the absence of strong Lyman alpha scattering of light at the appropriate redshifted 
wavelengths placed strong constraints on the amount of neutral hydrogen in the intergalactic 
medium (Shklovskij 1964; Scheuer 1965; Gunn & Peterson 1965). The absence of significant 
neutral opacity in a uniformly distributed component of the IGM is generally referred to as 
the Gunn- Peterson effect. 

While it remains conceptually possible that there is a diffuse substrate of the IGM in 
a multi-phase gas, modern interpretations of the Lya forest of neutral hydrogen absorption 
lines suggest instead that the IGM is a highly structured warm ionized gas in which the 
amplitude of the perturbations grows gravitationally. In this interpretation the Lya forest 
is produced by scattering from the very small neutral hydrogen fraction of this undulating 
density intergalactic gas (e.g., Cen et al. 1994; Zhang, Anninos, & Norman 1995). Regions 
of minimum H I optical depth, in this scenario, occur in cosmic mini-voids, which are un- 
derdense expanding regions where growth is essentially frozen (Meiksin 1994; Reisenegger & 
Miralda-Escudc 1995). The cosmic mini- voids are the closest analogs in this model to the 
Gunn-Peterson effect in a homogeneous diffuse gas, and the simple physics of these regions 
makes them powerful probes of the metagalactic ionizing flux relative to the baryon density 
of the universe. 

Because of the rapid increase in the gas density with redshift, it is expected that the 
fraction of the quasar spectrum returning to the continuum level will diminish rapidly at 
the higher redshifts. As a specific example, Zhang et al. (1997) find, at ^ = 5, for Hq = 
50 km s~^ Mpc~^, flf, — 0.06 and a Haardt-Madau (1996) spectrum and ionizing flux, that 
less than a fraction of a percent of the spectrum will have t(H I) < 0.1. Moreover, they find 
that the spectrum of a 2; = 5 quasar will have an average flux in the region between 1050 A 
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and 1170 A that is 10% of the continuum value which would be present in the absence of 
Lyman alpha scattering. 

The newly discovered z = 5 quasar J033829.31 + 002156.3 (Fan et al. 1999) presents 
an opportunity to look at this question at the highest redshift yet available, and we report 
here low and high resolution spectral observations of this object. As Schneider, Schmidt, 
& Gunn (1991b) found in the z = 4.897 quasar PC1247+3406, 1 - Da (the Oke index Da 
is defined in §3) is higher than the model values (0.25 in J033829.31 + 002156.3 and 0.36 
in PC1247+3406) and we also find that there are significant regions of low optical depth 
(r < 0.1). Our results are similar to the minimum optical depth of r = 0.02 ± 0.03 at 
z = 4.3 found in the z = 4.7 quasar BR1202-0725 by Giallongo et al. (1994) and to the 2 a 
lower limit of 0.1 found by Williger et al. (1994) in the z = 4.5 quasar BR1033-0327. Both 
the Oke index and the low minimum opacity require a relatively high ionizing flux at these 
redshifts. 

2. Observations 

Long-slit spectroscopic observations of J033829.31 + 002156.3 were obtained with LRIS 
(Oke et al. 1995) on UT 10 Feb 1999 with a 1^0 wide slit and 400£/mm grating blazed at 
8500 A (spectral resolution 8.1 A), and on UT 17 Feb 1999 with a 1'.'5 wide sht, with both 
the 400£/mm grating blazed at 8500 A (spectral resolution 12.3 A), and with the 300^/mm 
grating blazed at 5000 A (spectral resolution 17.3 A). Seeing was 0^'7-0^'9 FWHM on the first 
night, and 0'.'7 — I'.'O on the second night; observations were made under dark, photometric 
conditions in each instance. A sequence of three exposures, shifted by 10" along the slit, 
was taken in each configuration with the position angle of the slit set to the parallactic 
angle for the middle exposure. This yielded net integrations of 3600 s (5542-9350 A; 8.1 A 
resolution; I'.'O sht), 3600 s (5799-9609 A; 12.4 A resolution; 1'.'5 slit), and 3000 s (3673-8708 
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A; 17.3 A resolution; I'.'S slit). A GG495 blocking filter was used to suppress second order 
blue light. The observations are summarized in Table [l|. Calibration stars (HZ4 and HZ44) 
were observed at the parallactic angle for each slit/grating configuration and used to flux 
calibrate the data. The calibration stars used for the LRIS and HIRES data are Category 
1 HST Standard Stars (Bohlin & Lindler 1992), which are part of an ongoing effort to 
develop a consistent set of fundamental flux standards for observations with the Hubble 
Space Telescope over the wavelength region from 1050-10000 A based on the Oke (1990) 
measurements taken on the Double Spectrograph of the Palomar 5-m telescope for optical 
measurements. Although both HZ4 and HZ44 have long-wavelength flux calibration data 
measured by a variety of authors (e.g., Oke 1974; Stone 1977; Massey & Gronwall 1990), 
for self-consistency we use the Oke (1974, 1990) measurements. The estimated galactic 
extinction is Eb-v = 0.10 in this direction (Schlegel, Finkbeiner, & Davis 1998). 

The HIRES (Vogt et al. 1994) observations were obtained using the red collimator and 
the Dl decker (I'.'l x 14" slit; resolution R = 36, 000) on the nights of UT 14 and 15 Feb 1999 
under conditions of variable seeing. The four 2400 s exposures were taken at P. A. = 100° 
(non-parallactic angle) using the atmospheric dispersion corrector to maintain the guide 
object for this observation within the guide camera field. The white dwarf standard star 
G191-B2B (Oke 1990) was used to correct the continuum fit, which was then flux-calibrated 
by comparing an appropriately smoothed spectrum with the lower resolution LRIS data. The 
combined LRIS spectra, smoothed to the lowest resolution data taken with the 300£/mm 
grating (resolution 17.4 A), is shown in Figure 0, and the HIRES spectrum in Figure ^. Prior 
to any extinction corrections, both continua are well fit by a flat f^, spectrum in the line-free 
portions of the red spectrum, and this is shown by the dashed line. 

The fluxed LRIS spectra are in excellent agreement with the photometric measurements 
of Fan et al. (1999), and with their estimated ABi4^^q continuum magnitude of 20.01 (35 
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/iJy). A strong Mg II doublet is visible superposed on the C IV emission (Figure 0). We 
measure a doublet ratio = 1.25 and a redshift z = 2.304 for this system. If the strong 
rest equivalent width (M/a, rest (2796) = 1.73 A) indicates a large galaxy along the line-of-sight 
(Bergeron & Boisse 1991; Steidel, Dickinson, & Persson 1994), the quasar may be amplified 
by lensing. 



3. Discussion 

The simplest measures of the opacity of the neutral hydrogen are the Oke and Korycanski 
(1982) indices. Following Schneider, Schmidt, & Gunn (1991a) and Schneider et al. (1991b) 
we measure the index 

/i, (observed) \ 
(continuum) / 

in the rest-frame wavelength range 1050 A to 1170 A, obtaining = 0.75 in both the 
lower resolution spectra. This can be compared with the values of 0.55 to 0.74 measured by 
Schneider et al. (1991b) and Kennefick, Djorgovski, & de Carvalho (1995) in quasars in the 
range z = 4.35 — 4.5 and with the values of 0.74 obtained by Schneider et al. for the z = 4.733 
quasar PC1158+4635 and 0.64 for the z = 4.897 quasar PC1247+3406. It is considerably 
smaller than the value of 0.9 predicted by Zhang et al. (1997) for models with the Haardt 
and Madau (1996) Ju evolution, with Vlh = 0.06 and Hq = 50 km s^^ Mpc^^, suggesting 
that these models are underestimating J^/Ql at high redshift. 

This can also be examined by looking at the regions of minimum opacity in the much 
higher resolution HIRES spectrum. As can be seen from Figure 0, portions of the HIRES 
spectrum return close to the extrapolated continuum. 14% of the spectrum has r < 0.1 in 
the observed wavelength range 6800 — 7000 A if we extrapolate the continuum with a 
power law. This fraction is not highly sensitive to the power law choice, and reduces only 
to 11% if we use a z/^ extrapolation. The regions of minimum opacity around 6790 A and 
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6948 A have extremely low optical depth. The mean optical depth between 6946 A and 6949 
A is -0.05 ± 0.05 for a extrapolation and 0.05 ± 0.05 for a extrapolation, where the 
errors are based only on the statistical noise in the region. In the subsequent discussion we 
will adopt what we beheve is a conservative limit of r < 0.1 at this redshift of 4.72. 

For a uniformly distributed IGM and qo = 0.5, 

where h^o is the Hubble constant in units of 50 km s"^ Mpc~^, J_22 is the metagalactic flux 
at the Lyman edge (J^ = 10~^^ J_22(i^/t'9i2) " ergs cm~^ s^^ Hz^^ sr^^ with a = 0.7 ) and 
T4 is the temperature of the gas in units of 10"^ K (Giallongo et al. 1994). For r < 0.1 we 
obtain 

hi, < 8.4 X 10-3 hl-^ J% Tl-^'\ 
Using Haardt & Madau's value of J_22 = 0.4 at z = 4.72 would give Qfe/iio < 0-005 h^^ Tf^^^ 

If, however, the minimum opacity arises in cosmic minivoids, this number must be 
corrected upward to allow for the underdensity in these regions. Simulations with current 
cosmological parameters suggest that the proflle of the density distribution at 2; = 5 will 
be centered on a mode of an underdensity of about 0.5, and will have very little volume 
indeed at underdensities of 0.2 or less (e.g., Zhang et al. 1998). The underdense regions are 
relatively cold, with temperatures of around 5000 K. If we conservatively assume that the 
minimum opacity portions of the current spectrum arise in regions which are underdense 
by a factor of five, and set T4 = 0.5, we would obtain ^hhl, of 0.019 for J_22 = 0.4. We 
can compare this estimate of the baryon density to that derived from measurements of 
primordial D/H. The lowest value of fife/iso (0.02) that might currently be possible would 
correspond to D/H = 2 x 10^^ (Songaila, Wampler, & Cowie 1997; Webb et al. 1997) 
whereas D/H = 3.3 x 10-^ or fib hj, = 0.077 ± 0.006, is obtained by Tytler's group (Buries 
& Tytler 1998). Even with the minimum ^bhl, = 0.02, we require J_22 > 0.4 at 2; = 4.72, 
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whereas the low D/ H value would require J_22 > 7, more than an order of magnitude higher 
than the Haardt-Madau estimate. 

We thank T. Bida, R. Goodrich, G. Wirth, J. Aycock, C. Sorenson, and W. Wack 
for their assistance in obtaining the observations and M. Ranch and A. Meiksin for helpful 
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AST 96-17216 and NASA grant GO-7266.01-96A from the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS 5-26555. R.G.M. thanks the 
Royal Society for support. 
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Table 1. Spectroscopic Observations of J033829.31+002156.3 







^exp 


A A Range'^ 


Slit 


PA^ 


Resolution'^ 




FWHM^ 


IJT Date 


A.M. 




(A) 


(arcsec) 


(°) 


(A) 


Grating'' 


(arcsec) 


LRIS 


















10 Feb 1999 


1.08 


3600 


5542-9350 


1.0 X 274 


28.0 


8.1 


400/8500 


0.72 


17 Feb 1999 


1.12 


3600 


5799-9609 


1.5 X 274 


50.0 


12.4 


400/8500 


0.71 


17 Feb 1999 


1.30 


3000 


3673-8708 


1.5 X 274 


62.0 


17.3 


300/5000 


1.00 


HIRES 


















14 Feb 1999 


1.15 


2400 




1.1 X 14 


100.0 




Red 


3.0 


15 Feb 1999 


1.10 


7200 




1.1 X 14 


100.0 




Red 


1.2 



^The usable wavelength range for the multi-order HIRES data is ~ 6800 - 8000A 
''Set close to estimated parallactic angle at mid-exposure in the sequence of LRIS observations 
= 36000 for HIRES 

'^Identified by ruled £/mm and blaze wavelength in A for LRIS; for HIRES this is the RED 

configuration with Dl decker 

^Measured FWHM of quasar profile on slit 
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Table 2. ;s = 2.304 Mg II Absorber 





A 


Wx^ rest 


ID 


(A) 


(A) 



Mg II 2796 9236.7 1.73 
Mg II 2803 9260.4 1.38 
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Fig. 1. — Both the 300£/mm and 400£/mm fluxed LRIS spectra are shown, in //Jy versus 
wavelength. The 300^/mm spectrum is shown over the wavelength range 4000 - 8200 A and 
the combined 400^/mm data from 5700 A to 9700 A. The spectra agree extremely well in 
the overlap region. The dashed and dotted lines show and i/^ continuum fits, normalized 
to the region 7700 — 8000 A, which lies redward of the atmospheric A band but blueward of 
the Si IV emission line. 

Fig. 2. — The flux-calibrated HIRES spectrum over the range 6800 - 7700 A. The spectrum 
has been smoothed to a resolution of 2 A. The dotted and dashed lines show the same 
continuum fits as in Figure 1. 
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